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* Abstract 

The aeromechanics 1 stability of a helicopter 
in ground resonance was analyzed, by incorporating 
five .different aerodynamic models in the coupled 
rotor/fuselage analysis. The sensitivity of the 
results to changes in aerodynamic modelling was 
carefully examined. The theoretical results were 
compared with experimental data and useful con- 
clusions are drawn regarding the role of aerodyna- 
mic modeling on this aeromechanical stability 
problem. The aerodynamic model which provided the 
best all around correlation with the experimental 
data was identified. 

Nomenclature 

a - lift curve slope 


X - total inflow, X » X^-HSX ; also inflow 

mode^ in figures only 

X l ,X lc’ X ls ~ inflow variables 

6X - unsteady wake induced perturbational 

inflow 

(o - modal frequency, imaginary part of s 

H - rotor R.P.M. 

<t> - body roll mode 

ip - azimuthal angle or nondimens ional time 

- fit 

p - density of air 


C 


C 


1 


Slx'Sly 

[L] 


- lift deficiency factor 

- coefficient in inflow equations, C, - 
0.5 or 1.0 

- Thrust coefficient 

- moment coefficients in roll and pitch 

- induced flow matrix 

- mass flow rate 


O - modal damping, real part of s 

o - solidity ratio 

0 - body pitch mode 

0^ - collective pitch setting of the blade 

C lc ,5l - cyclic lag coordinates 

^P*^R ~ progressing and regressing lag modes 

respectively 


M - rotor aerodynamic moment 


1. Introduction 


Mj - apparent inertia 

r - radial location of a typical blade 

section 

R - rotor radius 

s - eigenvalue 

dT - differential thrust 

B lc ,B ls - cyclic flap coordinates 

0p,(3p - progressing and regressing flap modes 

respectively 

y - Lock number 

y* - equivalent, reduced or effective Lock 

number 

Xq - steady or mean inflow 

* 
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Unsteady aerodynamics have a significant in- 
fluence on the aeroelastic and the aeromechanical 
stability characteristics of helicopters. The 
mathematical sophistication of refined unsteady 
aerodynamic models is sometimes prohibitive to in- 
corporate in the aeroelastic analyses and therefore 
it is quite frequent that rotary-winp aeroelastic 
analyses are based upon quaFisteady aerodynamic 
theory. Fortunately, there are some relatively 
simple unsteady aerodynamic models, known as inflow 
models, which can be conveniently incorporated in 
the aeroelastic and aeromechanical studies of heli- 
copters. These simple models are based upon the 
definition of certain inflow parameters which re- 
present essentially the unsteady wake-induced flow 
through the rotor disk, A number of such inflow 
models are available in the literature; however the 
applicability of a particular model to a given rotor 
dynamic problem and the sensitivity of the stability 
boundaries to the choice of the inflow model and 
comparisons with experimental data have not been 
considered in detail in the literature. Bousman 1 
has carried out an experimental study of the aero- 
mechanical stability of a hingeless rotor supported 
on a special gimbal which simulated the pitch and 
roll degrees of freedom. The availability of this 
high quality experimental data provides an 
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opportunity; (a) to test the validity of mathemati- 
cal models representing the coupled rotor/fuselage 
dynamics and (b) co determine the influence of 
various aerodynamic models on this aeromechanicai 
problem. Bouaman attributed some of the discrepan- 
cies, found between the theoretical results present- 
ed in his paper and the experimental results, to 
dynamic inflow. This conclusion was examined by 
Johnson^, in a recent study, where the unsteady 
aerodynamic effects on the rotor were represented 
by a dynamic inflow model^. Johnson* results 
with the dynamic inflow model 3 * ^ indicated better 
agreement with the experimental data than the 
results obtained using the quasi-steady aerodynamic 
model. Using the coupled rotor/body model^»6 with 
simple quasi-steady aerodynamics, the authors^ 
also obtained good agreement with the experimental 
results generated by Bousman*. Based on the agree- 
ment with the experimental data, they concluded 
that the coupled rotor/fuselage model developed, 
was reliable. 


The purpose of this study is to extend Ref. 7 
and study the sensitivity of the results obtained 
to changes in the aerodynamic assumptions used. 

To accomplish this objective, five different aero- 
dynamic models were incorporated in the mathemati- 
cal model representing the coupled rotor/fuselage 
dynamics and the sensitivity of the stability 
boundaries to changes in aerodynamic modelling was 
determined „ The theoretical results were compared 
with the experimental data and based on this com- 
parison, conclusions are drawn regarding the 
selection of the parameters used in defining these 
aerodynamic models. 

2. Aerodynamic Models Used in the Analysis 

The aerodynamic models, incorporated in this 
aeromechanicai stability study representing a 
coupled rotor/fuselage system, were: (a) quasi- 
steady aerodynamics, (b) two different perturbation 
inflow models and (c) two different dynamic inflow 
models. A brief description of these aerodynamic 
models is provided below. 

Quasi-steady Aerodynamic Model 


along the rotor blade, in the blade fixed, rotating 
coordinate system. On the other hand, the inflow 
models represent the global effects of the un- 
steady wake and therefore they are applicable to 
the complete rotor. The various inflow models are 
described below. 

Perturbation Inflow Model 


Prior to describing the perturbation Inflow 
model, it is useful to clarify certain aspects of 
the terminology used in the literature which deals 
with this subject. In some cases, the perturbation 
inflow model is referred to as quasi-static inflow 
model^ and in other cases as quasi-steady inflow 
mod el 

The induced flow-field acting on a helicopter 
rotor affects both rotor equilibrium (trim load- 
ings) and rotor response (transient loading). 
Therefore, it is reasonable to assume that the 
induced flow will also be affected by the oscilla- 
tions of the rotor. This assumption is the basis 
of both the perturbation inflow models and dynamic 
Inflow models. A detailed derivation of these 
inflow models can be found in Refs. 3 and 9. 

In these models the total induced velocity on 
the rotor disk due to the wake is assumed to con- 
sist of two parts: (1) a steady inflow, (for 

trim loadings) and (2) a perturbation inflow, 51, 
(for transient loadings). Therefore, the total 
induced velocity normal to the rotor disk is ex- 
pressed as 


A « X Q + 6A (1) 

Assuming that the perturbation inflow, 6A, varies 
azirnuthally as well as linearly along the radius, 
the total inflow can be written as 

X - X Q + X 1 + X lc | cos* + X ls ^ sin* (2) 

where the inflow* variables X it X^ G , X^ g are func- 
tions of time. These inflow variables are related 
to the perturbational thrust, roll and pitch 
moment coefficients through the following relation. 


The quasi-steady aerodynamic model, employed 
in the analysis, is based on Greenberg’s® formula- 
tion of unsteady aerodynamic loads on an oscillatory 
airfoil in a pulsating flow, Greenberg* s theory 
Is a modified form of Theodorsen's unsteady aero- 
dynamic theory. The quasi-steady model is obtained 
by assuming C(k) - 1 and neglecting the apparent 
mass terms (noncirculatory terms). In this model, 
the assumption of C(k) • 1 implies that t unsteady 
wake effects are totally neglected 

In flow Models 

The inflow models represent the unsteady wake 
effects in a simple form. In these models, the 
unsteady wake-induced flow through the rotor disk 
Is defined by a set of inflow variables and these 
variables essentially provide a corvt :tion to the 
Inflow assumed in the quasi-steady aerodynami 
theory. When Inflow models are used in the analysis 
of rotor dynamic problems the blade loads have to 
be calculated from the quasi-steady aerodynamic 
expressions. An important fact to be noted is that 
the quasi-steady aerodynamic model is a two dimen- 
sional local model and hence it is applied at a 
typical cross section located at a spanwlse station 



where P.A stands for perturbational aerodynamics. 
The elements of [L] can be obtained either theoreti- 
cally, by using momentum theory 3 *^ or experi- 
mentally™. 

In ground resonance type of aeromechanicai 
problems, the inflow variable X^ does not couple 
with the body and cyclic blade degrees of freedom 
and hence it does not have to be considered in the 
analysis. Thus only the equations for the inflow 
variables X^ c and X^ s are relevant to this specific 
problem and these can be written as 
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For axial flow through the rotor, which cor- 
respondato the present case, the elements of fL] 
can be obtained by applying momentum theory?, I e 
differential thrust on an elemental area dA 
(- rdrdip) of the disk is related to the inflow by 
the equation 

dT » m 2XfiR (5) 

It should be mentioned that by relating the total 
differential thrust (steady and perturbation) to 
the total induced velocity (steady and perturba- 
tion) in the form, given in Eq, (5), it is assumed 
that the thrust-inflow relation is the same for 
steady as well as perturbational conditions. This 
basic assumption implies that the variation of the 
forces on the rotor is sufficiently slow so that the 
classical actuator disk theory is valid for both 
steady and perturbational inflow velocities. 
Therefore, this inflow theory is also recognized to 
be a low frequency approximation to the unsteady 
aerodynamics of the rotor. 

Following Johnson^, the mass flow rate in 
Eq. (5) can be written as 

m = pX^RdA (6) 

It is important to note that the mass flow rate 
is defined with respect to the steady or mean value 
of the inflow Aq. 

The aerodynamic pitch and roll moments on the 
rotor disk, acting at the hub, can be obtained by 
taking moments of the elemental thrust about the 
hub center and integrating over the complete rotor 
disk. The pitch and roll moments are 


where the mass flow is defined with respect to the 
total induced velocity X, the inflow equations for 
X^ c and X^ 8 become 


(13) 

P.A 

Comparing Eqs. (11) and (13), it is evident that 
depending on the definition of mass flow rate, i.e. 
Eq. (6) or (12), the coefficients of the elements 
of [L]“l matrix differ by a factor of two. 



Equations (11) or (13) are complete only after 
identifying the right hand side. This is done by 
obtaining expressions for the moment coefficients 
using blade element theory. Once these have been 
obtained, a relation is established between the 
Inflow variables and rotor blade motion. It was 
shown in Refs. 3, 9 and 11 that incorporation of 
a perturbational inflow model, as represented by 
Eq. (11) or (13), in rotor dynamic problems yields 
a modification of the aerodynamic loads acting 
on the blade which can be represented by a reduced 
or effective Lock number 

Y * = Cy (14) 

If Eq. (11) is used in the rotor dynamic problem, 
the lift deficiency factor C becomes 

C i (15) 


J *R r* 2 tt 

0 J 0 - r cosWT (7) 

/•R /*2tt 

“roll ■ J 0 J 0 r Si "* dT (8) 

Substituting Eqs. (2), (5) and (6) in Eqs. (7) and 
(8) and integrating , the pitch and roll moments 
become 


This factor is found to be equal to the low fre- 
quency approximation of Loewy’s lift deficiency 
function for harmonic loadings^. On the other 
hand, Eq. (13) produces a lift deficiency factor 


C 


1 


1 + 



(16) 


which is higher than that given in Eq. (15). 


pitch 

7T 

" 2 

P R 3 x 0 x lc (nR ) 2 

(9) 

*roll 

7T 

“ 2 

pR 3 X 0 X ls (fiR) 2 

(10) 


After nondimensionalization, the relation between 
inflow variables and the perturbational aerodynamic 
moment coefficients becomes 


The two perturbation inflow models, used in 
the present analysis, can be written in a general 
form as 




On the other hand, if the mass flow rate m, is de- 
fined as (following Peters and Gaonkar?) 

m - pXftRdA (12) 


when Ci - O.S.Eq. (17) corresponds to Eq, (11) and 
when Ci - 1.0 it corresponds to Eq. (13), 

The concept of equivalent Lock number in the 
coupled rotor/fuselage type problems appears to 
involve a certain inconsistency. The fuselage 
equations of motion in pitch and roll contain terms 
due to both aerodynamic hub moments and aerodynamic 
hub forces. When using the perturbation inflow, 
one can make the observation that only the Lock 
number associated with the aerolynamic moment terms 
is modified, however the Lock number associated 
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with the aerodynamic force terms remains unchanged. 
Of course, the reason for this inconsistency lies 
in the formulation of the equations for the inflow 
variables which are related only to the hub moments, 
as given in Eq. (4), 

Dynamic Inflow Models 

The perturbation inflow model does not account 
for the time lag between the aerodynamic load and 
the time variation in inflow* The dynamic inflow 
models represent an extension of the perturbation 
inflow model by taking into account the time lag 
between the aerodynamic loading and the response. 
When using the dynamic inflow model the equations 
for Xj c and Xj s can be written as 


model 2 

Case (d) : perturbation inflow model with C^*1.0 

Case (e): dynamic inflow model with Cj-1.0 and 

M 1 =0. 1 132 

These aerodynamic models can also be viewed 
as a special case of a general dynamic inflow 
model. When M^-0, the general inflow model be- 
comes a perturbation inflow mode] and when Mi + ®°, 
the effects of inflow perturbations are totally 
eliminated and the resulting model is a quasi-t^eacfy 
aerodynamic model. 

3, A Brief Summary of the Experiment 



where Mi represents the nondimensional apparent 
inertia associated with the inflow and the quantity 
Ci is either 0,5 or 1,0, depending on the definition 
of mass flow rate. The value of Mi can be obtained 
either theoretically or experimentally. Tuckerman 12 
evaluated the apparent inertia associated with an 
impermeable disk subject to an angular acceleration. 
The nondimensional value of the apparent inertia 
was found 12 to be Mi * 0.1132. This theoretical 
value is also supported by parameter identification 
studies 1 ^. In Ref. 13, it was noted that Mi can also 
be influenced by the pressure distribution on the 
rotor and hence Mi could be also assumed to be a 
function of rotor loading distribution. In Ref. 

14, the identified value of Mi is found to vary 
between 0,05 to 0,2. In the present analysis, the 
value of Mi is assumed to be the theoretically 
evaluated value i,e., Mi ■ 0.1132. The implication 
of using Eq. (18) in rotor dynamic problems 1 ^, 
under harmonic loadings, can be shown to be equi- 
valent to a modification of the Lock number, which 
can be written as 


Y* " Y 


1 - 


16C-I \ n 16M n (i<o) 

H — + 


(19) 


oa 


oa 


* yC 


Equation (19) Indicates that addition of an apparent 
Inertia term to the perturbation Inflow model intro- 
duces a phase lAg between the aerodynamic loads and 
the response. Furthermore the value of C Is now 
different from the previous values given in Eq, (15) 
and (16). 


The five aerodynamic models, described briefly 
above, were selected for Incorporation In this 
study. Using these theories, the sensitivity of 
the aeromechanlcal stability problem to changes in 
the aerodynamic assumptions was Investigated. For 
convenience, these five aerodynamic models are 
concisely summarized below: 

Case (a): quasl-steedy aerodynamics 

Case (b): perturbation Inflow model with Ci-0,5 

Case (c): dynamic Inflow model with Cj-0,5 and 

Mi«0.1132 which corresponds to Johnson's 


A clear description of the experimental set 
up, used for simulating the fundamental aspects of 
the aeromechanlcal stability of a hingeless rotor 
helicopter, was presented in Ref. 1. The rotor 
consisted of three blades and five different con- 
figurations were tested. The different configura- 
tions represent different blade parameters char- 
acterized by the nonrotating natural frequencies 
of the blade in flap and lag, pitch-lag coupling 
and flap-lag coupling. The rotor was designed 
such that most of the blade flexibility is concen- 
trated at the root by building in root flexures. 

The rotor assembly was supported on a gimbal 
which had pitch and roll degrees of freedom. In 
this paper the analytical results obtained are com- 
pared with the experimental results, presented by 
Bousman, for rotor configurations 1 and 4, where 
the designation of these configurations is consist- 
ent with those in Bousman' s paper 1 . A brief 
description of these configurations is presented 
for the sake of completeness, additional informa- 
tion can be found in Refs. 1, 2 and 7. Configura- 
tion 1 had different stiffnesses in flap and lag 
respectively; the corresponding nonrotating flap 
frequency was 3.13 Hz and that for lead-lag was 
6. /0 Hz. Configuration 4 was a matched stiffness 
case where the nonrotating flap frequency was 6,63 
Hz and that for lead-lag was 6,73 Hz. The pitch- 
flap and pitch-lag coupling for these two config- 
urations was zero. For cases where the pitch 
angle was nonzero, the experimental rotor was 
designed such that pitch changes were introduced 
outboard of the flexures and therefore the struc- 
tural flap-lag coupling for these cases was zero. 
The blade was also designed to be very stiff in 
torsion. 


4 , Method of Solution 

The degrees of freedom considered in this 
aeromechanical stability analysis are: the funda- 

mental flap and lag modes of the blade and the 
pitch and roll degrees of freedom of the body. In 
this class of problems, it has been established 
that the collective flap and lag modes do not 
couple with the body motion and thus, these modes 
are not considered. Since the inflow variable Xj 
also has the role of a collective mode, it need 
not be considered, Therefore, the total number of 
degrees of freedom governing the aeromechanical 
problem are six. They *re: cyclic flap (fli c ,B]. g ), 

cyclic lead-lag (Ci^,4i s ). body pitch (0) and 
body roll ( 40 . For the cases when the dynamic 
Inflow models are used, two additional degrees of 
freedom, namely Xx c and X j_ 8 , are also present in 
the problem. 
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The solution of the coupled rotor/fuselage 
problem follows essentially the procedure outlined 
in Refs. 6 and 7. The procedure for obtaining the 
stability boundaries of the system consists of the 
following steps. 

1 . Evaluation of the equilibrium position of the blade. 

2. Linearization of the nonlinear equations of 
motion about the equilibrium position. 

3. Transformation of the linearized equations with 
periodic coefficients to equations with constant 
coefficients by using a multiblade coordinate 
transformation. 

4. Evaluation of the eigenvalues of the linearized 
system with constant coefficients to obtain the 
stability boundaries. 

The eigenvalues appear in complex conjugate 
pairs, s * a+ la). The real part of the eigenvalue 
represents the modal damping and the imaginary part 
modal frequency, respectively. The mode is stable 
if 0 is negative and it is unstable if a is posi- 
tive. 

In the present problem, the number of complex 
eigenvalue pairs depends on the type of aerodynamic 
model used in the analysis. When quasi-steady 
aerodynamics or the perturbation inflow models are 
used, there are only six pairs of complex eigen- 
values, each one representing one of the six degrees 
of freedom, namely, Bi c , B ls , C lc> 4 ls , 8 and 
The modes corresponding to the rotor degrees of 
freedom (B ic , 8x s , Ci c , Cj, s ) are referred to either 
progressing or regressing mode depending on the 
numerical value of the rotating natural frequency. 

A more detailed description of this terminology 
can be found in Refs. 3 and 7. When the dynamic 
inflow model is used, the six eigenvalue pairs are 
augumented by one additional pair of eigenvalues 
corresponding to the inflow variables. Since the 
equations for the inflow variables Ax c and Ai s are 
given in first order state variable form, Eq. (18) 
the stability analysis will yield only one pair of 
eigenvalue corresponding to these two inflow vari- 
ables, The mode corresponding to this eigenvalue 
pair is designated as the "inflow mode", (A), 
following Johnson’s 2 terminology. 

5, Results and Discussion 

In the present study, aimed at predicting the 
aeromechanical stability of a model helicopter, the 
behavior of the model is studied at various values 
of rotor speed ft. Two rotor configurations are 
analyzed. Configuration 1, in which the nonrotating 
flap frequency is lower than the nonrotating lag 
frequency, and configuration 4, in which these two 
frequencies are almost equal, which corresponds to 
a matched stiffness configuration. These different 
configurations have an influence on the dynamic 
behavior of the coupled rotor/fuselage system. In 
a matched stiffness configuration the structural 
flap- lag coupling is eliminated. Furthermore the 
root torsional moment due to the combined flap-lag 
motion, which is somewhat similar to an effective 
flap-pitch and lag-pitch coupling, is also elimin- 
ated. It should be mentioned however that these 
effective flap-pitch and lag-pitch couplings are not 
structural couplings. It was mentioned in the 
previous section that the experimental model was 
designed so as to eliminate structural flap-lag 
coupling, for these configurations. Therefore, the 
difference between these two configurations con- 
sists of the root torsional moment due to combined 
flap-lag motion which is present in configuration 1 


and absent in configuration 4, This root moment 
acts as an exciting moment for the body pitch and 
roll motions. 

The numerical data used in the analysis is 
presented in the Appendix. It should Be mentioned 
that the roll inertia used in the present calcu- 
lations is slightly higher than the value 
(183 gm*m 2 ) provided in Ref. 1. The value for roll 
inertia used in our calculation is 104 gm«m 2 which 
is 6% higher than 183 gm*m 2 . This value of roll 
inertia was obtained by using the body spring 
stiffness in roll, provided by Bousman-^, such that 
the calculated nonrotating coupled roll frequency 
matches the measured frequency. 

5.1 Results for Configuration 1 

The results for Configuration 1 are presented 
in Figs. 1-8, The variation of various modal fre- 
quencies with ft are shown in Figs, 1-2, together 
with the experimental data, taken from Ref. 1. It 
can be seen from Fig. 1 that the analysis witfc 
quasi-steady aerodynamics predicts the modal fre- 
quencies which are in excellent agreement with 
the experimental results. Figure 2 presents the 
calculated modal frequencies for Case (b) , pertur- 
bation inflow with Ci»0.5, and Case (c) dynamic 
inflow model with Ci-0.5 and M],*0.1132, With the 
perturbation inflow model, the predicted frequen- 
cies for roll (4>) and pitch (9) are over estimated 
in the range ft > 300 R.P.M.. On the otherhand, 
by incorporating a time delay in the inflow model. 
Case (t), the calculated pitch and roll frequencies 
are in good agreement with the measured values. 
However, the predicted pitch frequency is still 
slightly higher in the range ft > 300 R.P.M. A 
similar trend was also observed in the results for 
Case (d) , perturbation Inflow model, with Ci*1.0, 
and Case (e) dynamic inflow model with Cx*1.0 and 
Mi-0.1132. 

It was mentioned earlier that the analysis 
with dynamic inflow model produces an additional 
eigenvalue corresponding to the inflow mode 0 ) . 

For Case (c) , there are two eigenvalues with 
frequencies below 0.6 Hz in the range ft > 200 R.P.M. 
as evident from Fig. 2. The frequency correspond- 
ing to one mode remains almost constant C«0.5 Hz), 
while the other decreases to zero and then increas- 
es. It is difficult to identify which one of 
these two corresponds to the flap regressing mode 
(8 r) and which one should be associated with the 
Inflow mode (A). The mode with the constant fre- 
quency, in Fig, 2, is identified as Inflow mode 
(A) and the other mode is identified as flap re- 
gressing mode (8r) . Johnson 2 also identified the 
mode with a constant frequency as inflow mode (A) 
and the second mode as flap regressing mode (8 r ). 
Some additional comments on this identification 
procedure will be made later. 

Figure 3 presents the variation of damping in 
flap regressing mode (Br) and Inflow mod t (A) with 
ft. It is evident from Fig. 3 that the damping in 
the flap regressing mode increases rapidly with 
ft for the analysis with quasi-steady aerodynamics. 
The introduction of the perturbation inflow model 
with Ci*0.5, Case (b) , drastically reduces the 
damping in 0 R mode. This reduction in damping is 
caused by reduced aerodynamic damping with pertur- 
bation inflow. For this case the relevant 
quantities are: solidity ratio 9 ■ 0,0494; lift 

curve slope a • 5,73 and steady inflow Xo*0,014, 
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Therefore the deficiency function C f based on Eq. 

(15), is C-C.284. Hence, the effective Lock number 
f*-0.284Y. This shows that perturbation inflow 
reduces the magnitude of aerodynamic forces by 
approximately 72%. In the case of dynamic Inflow 
with Ci-0,5 and Mi-0.1132, Case (c), the damping 
in the mode which is identified as the inflow mode 
(X) remains relatively low, but the damping in 
flap regressing mode (ftp) Increases with ft. These 
results indicate that the damping in flap regress- 
ing mode reverts to the value obtained in the 
analysis with quasi-steady aerodynamics, as a con- 
sequence of the time delay present in dynamic inflow 
model. This seems to contradict the earlier results 
published in Refs. 9, 10 and 17. It was mentioned 
in Ref. 17 that flap regressing mode damping is 
substantially decreased by dynamic inflow for small 
values of collective pitch setting of the blade, 
Furthermore, it was found in Ref. 10 that dynamic 
inflow reduces the damping in flap regressing mode. 
This raises e question whether the Inflow mode 
identified in Fig. 2, and also identified as such 
by Johnson?, is a flap regressing mode and the mode 
identified as the flap regressing mode is really 
an inflow mode. To ascertain the reliability of 
this identification procedure, the eigenvectors 
corresponding to these modes were also analyzed. 

Table I shows the eigenvectors corresponding to 
the mode identified as the flap regressing and 

the inflow mode (X) at ft - 900 R.P.M. It can be 
seen that in the flap regressing mode, the flap 
motion has a higher participation factor than the 
inflow variables. In the inflow mode, the flap and 
inflow variables have almost equal participation 
factor. Also in this (X) mode, the pitch and roll 
motions have substantial participation factors. 
However, from these results, one can conclude that 
the flap regressing mode and Inflow mode are highly 
coupled modes. 

Figures 4 and 5 illustrate the variations of 
damping in pitch as a function of ft. Using quasi- 
steady aerodynamics, a higher damping, in the range 
200 < ft < 800 R.P.M. is predicted as evident from 
Fig. 4. However when dynamic inflow. Case (c^, 
with C^-0,5 and Mi«0.1132 is used the damping is 
predicted very well in the range 200 < ft < 800 
R«P*M.> however the damping is somewhat under pre- 
dicted beyond ft ■ 800 R.P.M. Figure 5 shows that 
using the dynamic inflow model, Case (e) with 
Ci*1.0 and Mi«0,l 1 32, the damping predictions are 
in very good agre nt with the experimental results 
over the complete range of ft. When the value of 
Ci is inc -eased from 0.5 to 1.0, in the dynamic in- 
flow models, one finds that the corresponding damp- 
ing in pitch increases by 10% — 25% for ft > 400 
R.P.M. It is also evident from Figs. 4-5 that 
perturbation inflow models do not seem to predict 
the correct damping Jevels. 

The variation of damping in roll mode is shown 
in Figs. 6-7, When using quasi-steady aerodynamics, 
the damping, in the range of ft > 500 R.P.M,, is 
over predicted as evident from Fig. 6, Using the 
perturbation inflow model, with Ci-0.5, Case (b) , 
the damping in the range of ft < 600 R.P.M. Is under 
estimated, beyond this range of ft, the predictions 
are good. The damping levels predicted using the 
dynamic inflow model, with Ci«0.5 and Mi-0.1132, 

Case (c), sre In good agreement with experimental 
data over the whole range of ft, as shown in Fig. 6, 

It can be seen from Fig. 7 that calculations with 
the perturbation inflow model, with Ci«1.0, Case 
(d) , yield damping levels which are too low in the 
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range ft < 500 R.P.M. and too high for the range 
ft > 600 R.P.M. Calculations with the dynamic in- 
flow model, with Cpl.O and Mi-0.1132, Case (e) , 
also overpredict the damping in the range ft > 700 
R.P.M. It should be mentioned, that in the range 
ft > 700 R.P.M., when using dynamic inflow models, 
the percentage increase in roll damping as a re- 
sult of increasing from 0.5 to 1.0 is 5% to 12%. 

Based on the results obtained for the damping 
in the pitch and roll modes, it appears that the 
theoretical results are quite sensitive to the 
value selected for C\. It is also evident that 
introduction of a time delay in the inflov 1 model 
seems to be an important factor. Based on the over- 
all agreement with the experimental data, it ap- 
pears that the dynamic inflow model with CpO.S 
and Ml*0.1132, Case (c) , seems to be somewhat 
superior to the dynamic inflow model with C^l.O 
and Mj“0.1132, Case (e) . 

Figure 8 presents the variation of regressing 
lag mode damping with ft. The predicted damping 
levels are in good agreement with the experimental 
results in the range ft < 700 R.P.M. and ft > 900 
R.P.M. for all the aerodynamic models used. For the 
cases analyzed with perturbation inflow models. 

Case (b) with Ci a 0.5 and Case (d) with Ci=1.0, the 
value of ft at which the resonant peak occurs is 
shifted from 760 R.P.M. to 800 R.P.M. This shift 
is associated with the fact that when using both 
models the roll frequency predicted is higher than 
the experimental result and as a consequerce, the 
resonance is also shifted to a higher value of ft. 
Calculations with quasi-steady aerodynamics predict 
correctly the value of ft at which resonance occurs, 
however the stability of this mode is over predicted. 
In the analyzes with dynamic inflow models, the 
predicted damping levels are in excellent agreement 
with the experimental results, including the damp- 
ing at resonance. The level of agreement with 
experimental data found in this case is somewhat 
better than those shown in Refs. 1 and 2. This 
result seems to support the statement made in Ref. 

7, that the coupled rotor/fuselage model derived 
in Ref. 5 is a reasonably accurate model for the 
configuration tested in Ref. 1. 

5,2 Results for Configuration 4 

The results for Configuration 4 are presented 
in Figs, 9-16. Figures 9 and 10 show the variation 
of modal frequencies as a function of ft. It can 
be seen from Fig. 9 that all the frequencies except 
the one corresponding to 0.7 Hz are predicted well 
by the quasi-steady aerodynamic model. In the 
range 250 <ft< 350 R.P.M., the pitch, roll and flap 
regressing modes undergo a change in their charact- 
eristics. The flap regressing mode (8r) becomes 
a roll mode ($) and roll mode (<J>) becomes a pitch 
mode (0) and the pitch mode (0) becomes a flap 
regressing mode (8r). In this range of ft, the pre- 
dicted roll frequencies are higher than the measur- 
ed values. Quasi-steady aerodynamics is incapable 
of predicting a frequency close to 0.7 Hz in the 
range 300 <ft< 1000 R.P.M. Figure 10 illustrates 
the results for the cases where the perturbation 
inflow model, Case (b) with Cj-0.5 and dynamic 
inflow model, Case (c) , with CpO.5 and Ml*0.1132, 
were used. Although calculations based on the 
perturbation inflow model are capable of predicting 
a frequency close to the experimentally measured 
frequency of 0.7 Hz, the pitch and roll modes 
frequencies are overpredicted. With dynamic inflow 
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model, all the frequencies are predicted well. 
However, in the range 250 < ft <350 R.P.M., the roll 
mode frequency is still overpredict d. In this 
range of ft, none of the aerodynamic models used in 
the present study, is capable of correctly predict- 
ing the roll frequency. Johnson 18 attributed this 
discrepancy to either the deficiency of the aero- 
dynamic model or the presence of some higher mode 
of the rotor or body . In any case, this problem 
remains unresolved. In the range ft > 400 R.P.M., the 
mode with a frequency close to 0.7 Hz is Identified 
as inflow mode (A) and the other mode with a fre- 
quency which is lower than 0.7 Hz is identified as 
flap regressing mode (@r) < This identification is 
based on the analysis of the eigenvectors correspond- 
ing to these two modes. Table II shows the eigen- 
vectors of the identified inflow mode and flap 
regressing mode, for ft * 900 R.P.M. In this case, 
as in Configuration 1, these two modes are highly 
coupled. However, in the Br mode, the flap 
metion has a higher participation factor than the 
participation of the inflow variables. In the A 
mode, the body pitch has the highest participation 
factor, with the flap, body roll and inflow having 
almost equal participation factors. As a result 
of this identification procedure one finds that 
using the dynamic inflow model the damping for the 
flap regressing mode is predicted to have a value 
comparable to that obtained when using quasi-steady 
aerodynamics. This seems to contradict some 
results which have been published earlier in Refs. 

9, 10 and 17 where it was found that using dynamic 
inflow yields a substantial reduction in regressing 
flap mode damping. 

The variation of roll damping as a function of 
ft is illustrated in Figs. 11-12. Calculations 
based on quasi-steady aerodynamics overpredict the 
damping in the range ft > 300 R.PM. , as evident in 
Fig. 11. Calculations based on the perturbation 
inflow model, with Ci»0.5, under predicts the damp- 
ing in the range ft < 800 R.P.M. T he damping levels 
predicted by using the dynamic inflow model, with 
Cp0.5 and M]*0.1132, are in reasonable agreement 
with the experimental values. From Fig. 12, it 
is evident that using the perturbation inflow model 
with Cj-1.0, yields a damping prediction which is 
too high for ft > 700 R.P.M. Using the dynamic in- 
flow model, with Cj-1.0 and Mj»0.1132, yields damp- 
ing level predictions which are in reasonable agree- 
ment with the experimental results. 

The variation of damping for the pitch mode 

(0) and the mode which has been identified as the 
inflow mode (A) are presented in Figs. 13 and 14, 

It is evident from Fig. 13 that predictions based 
on quasi-steady aerodynamics yield higher values 
of damping than the measured valuss. Calculations 
based on dynamic inflow, with Cj-0.5 and Mi*0.ll32, 
predict the pitch damping well, but the doping 
associated with the inflow mode (A) is lower than 
the experimental values. Figure 14 shows t'iat 
using the dynamic Inflow model, with Cj-1.0 and 
M^O.1132, yields a higher value of pitch damping 
than measured in the test. The damping in inflow 
mode is also higher than the experimental values. 
From the results shown in Figs. 13 and 14, it is 
evident that an Increase in the value of C| from 
0.5 to 1,0 Increases the pitch mode damping by 
about 25%. Therefore, it can be concluded that for 
Configuration 4, as well as Configuration 1, the 
predicted damping level? in pitch and roll modes 
are quite sensitive to the dynamic inflow model 
used in the analysis. By using a different 


combination of the values of Cj and Mj, bettsr cor- 
relation wlt> tbe experimental results could have 
been nchleven. However, among the aerodyramic 
models emp" vyed in the present study, the dynamic 
inflow with Ci"0. c and M^*0.1132, Case (c) , 

yieli c _.er agreement with tne experimental 
result than the other aerodynamic models, 

Thr. variation of lag regressing mode damping 
with ft: is shown in Fig. 15. The resonant frequency 
obtained >"h tue perturbation inflow models ex- 
hibits a sir "t t ) a higher value of ft than the one 
observed in the experiment. Calculations based 
quasi-steady aozodviamics predict the damping 
levels and the resonant f r rque*tcy very we J l, * 
the mode is more stable at ivso r ance, than te- 
stability indicated In the lest. Calculation - 
with the dynamic inflow models yield vaults * ni'-i* 
are in excellent agreement with the experimtn. i 
data. 

Figure 16 shows the variation of regressing ’ag 
mode damping as a function of collective pitch 
setting of the blade, at ft = 1000 R.P.M. Calcu- 
lations with the dynamic inflow model, using 
C]«0.5 and Mi*0.1132, yields damping values which 
are in very good agreement with the measured values. 

6 , Conclusions 

The aeromechanical stability of a helicopter 
in ground resonance is analyzed, using five differ- 
ent aerodynamic models and the analytical results 
are compared with the experimental results. Based 
on the comparison, the following conclusions can 
be drawn: 

(1) For the aeromechanical stability problem 
studied here, .he perturbation Inflow models do 
not predict correctly the modal frequencies and 
damping. 

(2) Quasi-steady aerodynamics ptedicts the modal 
frequencies very well for Configuration 1, but is 
incapable of predicting a frequency of 0.7 Hz 
measured in the experiment, for Configuration 4. 

The damping in body roll and pitch modes are over 
predicted. The regressing lag mode damping Is 
predicted well. 

(3) The dyiamic inflow models predict the modal 
frequencies and damping values which are in very 
good agreement with the experimental results. 

This implies that for the present problem the time 
lag is an important ingredient in the dynamic in- 
flow model. 

(4) The predicted damping levels for the lag re- 
gressing mode, using dynamic inflow models, are in 
excellent agreement with the experimental results 
including the value at resonance. Tnis indicates 
that the mathematical model for the coupled rotor' 
fuselage system is accurate. 

(5) From the cases studied with dynamic inflow 
models, it is evident that the predicted damping 
levels for the body modes, increase when Cy is 
increased from 0.5, Case (c) , to 1.0, Case (e). 

For both rotor configurations, the pitch damping 
increases by 10% to 25% and the roll damping in- 
creases by 5% to 12%. 

(6) Based on the comparison of results obtained 
with various aerodynamic models, it seems that the 
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dynamic inflow model with C^-0.5 and M^-0.1132 is 
the most suitable aerodynamic model, for the speci- 
fic aeromechanical problem studied in this paper* 

(7) Identification of the flap regressing mode 
and the inflow mode proved itself to be quiLe 
compiler ed* These modes were identified by using 
the frequency information together wi h a careful 
examination of the eigenvectors. The results based 
on this identification procedure seem to indicate 
that when u^ing the dynamic inflow model the pre- 
dicted values cf damping for the regressing flap 
mode are comparable to those obtained with quasi* 
steady aerodynamics. 
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Number of blades 3 

Radius, cm 81.1 

Chord, cm 4.19 

Hinge offset, cm 8.51 

Blade Airfoil NACA 23012 

Profile drag coefficient 0.0079 

Lock number 7.37 

Solidity ratio 0,0494 

Lift cu*ve slope 5.73 

Height of rotor hub above glint) il, cm 24,1 

Blade Mans Properties 

Blade mass (to flap flexures), gm 209 

Bl.^de mass centroid (Ref. flexure 

centerline), cm 18.6 

Blade flap inertia (Ref. flexure 

centerline), gm«ra^ 17.3 


Blade Frequency and Damping 

Coni. 1 Conf. 4 


Nonrotating Flap frequency, Hz 3 . 13 6.63 

Nonrot&tlng Lag frequency, Hz 6,70 6,73 

Damping in lead-lag (2 critical) 0.522 0.532 

Body Mass Properties 

Rotary inertia in pitch, gm*m 2 633 

Rotary inertia in roll, gm.m 2 194 

Body Stiffness and Damping 

Pitch stiffness, gm*cm2/sec 2 «rad 0,8687x10® 

Roll stiffness, gm*cin /secured 0.1113xl0 10 

Danping in roll (2 critics. 1 ) 0.9292 

Damping in pitch (2 critical) 3.202 

Pitch frequency, Hr 1,86 

Roll frequency, Hz 3.81 
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— • — QUASISTEADY AERODYNAMICS 

PERTURBATION inflow, C, - 0.5 

37^} DYNAMIC INFLOW. Cy • 0.5. M-, - 0.1132 



Fig. 3 Variation of Damping in Regressing Flap 
Mode and Inflow Mode with ft;8 c «0; 
Configuration 1 


□ EXPERIMENT (Rtf. 1) 

— •— PERTURBATION INFLOW. Cy - 0.5 

OUASISTEADY AERODYNAMICS 

DYNAMIC INFLOW, C, - 0.5, M, - 0.1132 



Fig. 4 Body Pitch Mode Damping as a Function of U; 
9 c -0; Configuration 1 


0 EXPERIMENT (Rtf. 1i 

PERTURBATION INFLOW, Cj - 1.0 

DYNAMIC INFLOW, C, - 1.0, M, • 0.1132 
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0 EXPERIMENT (Rtf. 1) 

— •— PERTURBATION INFLOW, C, - 0 5 

QUASISTEADY AERODYNAMICS 

OYNAMIC INFLOW, C 1 * 0.5. M, - 0 1132 



Fig. 6 Body Roll Mode Damping as a Function of Q; 
0 c -O; Configuration 1 


DYNAMIC INFLOW, C, - 1.0, My - 0.1132 

PERTURBATION INFLOW. C, « 1.0 



Fig. 7 Body Roll Mode Damping as a Function of ft; 
G c »0; Configuration 1 
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Fig. 5 Body Pitch Mode Damping as a Function of fi; 
0 c **O; Configuration 1 
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o EXPERIMENT (Rtf. 1) 

— •— PERTURBATION INFLOW, C 1 - 0.5 

QUASISTEADY AERODYNAMICS 

DYNAMIC INFLOW, C, ■ 0.5. M, - 0.1132 
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O EXPERIMENT (Rtf. 1) 

DYNAMIC INFLOW. C, - 1.0. M, - 0.1132 

PERTURBATION INFLOW, C, - 1.0 


Fig. 10 Variation of Modal Frequencies with Cl; 
0-0; Configuration 4 
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Fig. 8 Regressing Lag Mode Damping as a Function 
of ft; 6 C “0; Configuration 1 
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Fig. 11 Variation of Roll Mode Damping with ft; 


6 c *0; Configuration 4 
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Fig. 9 Variation of Modal Frequencies with ft; 
0 c «O; Configuration 4 
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Fig. 12 Variation of Roll Mode Damping wifh ft; 

0 -0; Configuration 4 
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□ A EXPERIMENT (R«f. 1) 

QUASISTEADY AERODYNAMICS (0) 

DYNAMIC INFLOW | 

— • — DYNAMIC INFLOW (A) j 


C, « 0.5. M, - 0.1132 



Fig. 13 Variation of Pitch (0) and Inflow ( \ ) 
Modes Damping with ft; 0 c *O; 
Configuration 4 
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Fig. 15 Regressing Lag Mode Damping as a Function 

of ft; 9 *0; Configuration 4 
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Fig. 16 Regressing Lag Mode Damping as a Function 
of 6 c ; Configuration 4 
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Fig. 14 Variation of Pitch (0) and Inflow (A) 
Modes Damping with ft; 0 c *O; 
Configuration 4 
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DISCUSSION 
Paper No. 14 

INFLUENCE OF VARIOUS UNSTEADY AERODYNAMIC MODELS ON THE AEROMECHANICAL STABILITY 
OF A HELICOPTER IN GROUND RESONANCE 
P. Friedmann 
and 

C. Venkatesan 


Bill Bousman. U.S. Army Aeromechanics Laboratory : Jack Landgrebe put a perspective on dynamic 

inflow as not really contributing to the loads problem and I think that the conclusions from the 
first two papers show that although it has a strong effect on the flapping degree of freedom, 
from fie designer's point of view, it is not really important for the lag degree of freedom. I 
guess my question is to all three of these guys. Are there applications for dynamic inflow in 
something like the handling qualities area where simulation needs the speed of the model and has 
any work been done in here or are there paths that we should be going? 

Fr. *dm« nn: 1*11 tell you I was expecting this question so I have a slide. "Can I have the 
slide, please?" The last slide [Fig. 3] is something which in your experiment you might have 
data, but it wasn't in your paper so I don't know whether you have data or not. It shows the 
flap regressing mode damping with various kinds of aerodynamics. What it really shows you is 
that the damping with quasisteady aerodynamics is here. If you put in the perturbation inflow 
It knocks down this damping in the flap mode very significantly. And uhen you put in the 
dynamic inflow with of 0.5 it brings it up again to almost where the damping was with the 
quasisteady aerodynamics. You can also see how the damping of inflow mode changes as a function 
of 0. So in relation to the first question which you have asked I think that the better test 
for how mucn global t~uth is in dynamic inflow should really be based on the behavior of the 
flapping uode as has teen indicated by both Dave Peters and [Gopal] Gaonkar. Maybe in the 
future some calculations associated with that type of exaounation could be revealing. 

Dave Peters. Washington University : On the question about handling qualities, I think it 

definitely has an effect. There was one figure in the paper we didn't show which shows the 
pitch and roll moment on the rotor due to a roll oscillation or a pitch oscillation. As you go 
to an in of zero the slope of that curve then is the roll rate or pitch rate moment as a 
function of e or e , like a control derivative. There's more than a factor of two 
difference with or without dynamic inflow; almost a factor of three in one case. I think if you 
are going to do handling qualities, anything in the less than once per frequency range then the 
dynamic Inflow It going to have importance. That's a great paper, Peretz. I think we snould 
have an altar <ull and everyone who wants to put dynamic inflow in their analysis should come 
forward or something after a paper like that because it's really good. One question I had — on a 
model like yours, how much extra complexity does it take to put the dynamic inflow in? Is it 2% 
or 10%? Maybe you can give us a feeling for that. 

Friedmann : It may be 10> additional work. It's not really very difficult to do. Particularly 

if you have somebody as good as Venkatesan who does it. 

Bob Loewy, Rensselaer Polytechni c Institute : My question pertains to the off-diagonal terms in 

the L matrix and really takes up a little bit on Euan Hooper's earlier question on the earlier 
paper. And that is it seems to me that swirl would make those kind of terms nonzero and that 
particularly in tilt rotors and high speed forward flight you would expect more swirl than we 
are used to. I wonder If you have thought about these things? 

Friedmann : I have to phrase this very carefully. I am essentially somebody who uses dynamic 

inflow. I am not a person who has ambitions of improving dynamic inflow. I am a believer in 
unsteady aerodynamics and as a consequence you might be aware a year ago one of my students 
completed an arbitrary motion type of unsteady airfoil theory in which you can essentially do 
the same things you do with dynamic inflow, but for hover and forward flight. It is based on 
essentially an assumed wake. [It has] all the mathematical complexities and maybe mathematical 
fundamentals which an unsteady aerodynamics theory provides you and you don't have to use the 
assumptions which are embedded in dynamic inflow and cannot be removed. We have used this 
particular arbitrary motion theory to essentially extend the so-called Loewy lift deficiency 
function, which you might be familiar with, to arbitrary motions. In that AIAA paper which was 
given last year we have not been very successful. But since then Dr. Venkatesan has managed to 
do an arbitrary motion approximation to the theory and that theory can probably be used to 
capture the same behavior which is predicted by the dynamic inflow model and you might be able 
to see whether based on such a theory you do get off-diagonal terms or not. 
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